* C:\Dropbox\Rehm\Data\CNEF\cnef_1970-2014\slid\slid_rehm_v9_log_for-stata.do

cd C:\Dropbox\Rehm\Data\CNEF\cnef_1970-2014\slid\

**************
* ESIs, all ages
**************
* tabstat *ESI*_nYe [aw=weight], by(year) f(%4.3f)

set more off

clear
preserve
	input year ESI25_nYe  ESI50_nYe  ESI10_nYe  ESI00_nYe  upESI25_nYe  upESI50_nYe  upESI10_nYe  upESI00_nYe  gESI25_nYe  gESI50_nYe  gESI10_nYe  gESI00_nYe  upgESI25_nYe  upgESI50_nYe  upgESI10_nYe  upgESI00_nYe  
		 1994 0.135 0.069 0.245 0.462 0.105 0.036 0.257 0.538 0.119 0.041 0.239 0.462 0.121 0.054 0.265 0.538
		 1995 0.119 0.052 0.240 0.501 0.120 0.047 0.262 0.499 0.104 0.032 0.234 0.501 0.135 0.065 0.269 0.499
		 1996 0.107 0.044 0.219 0.457 0.113 0.041 0.260 0.543 0.090 0.028 0.213 0.457 0.128 0.057 0.268 0.543
		 1997 0.112 0.048 0.213 0.454 0.118 0.044 0.274 0.546 0.097 0.030 0.207 0.454 0.136 0.064 0.283 0.546
		 1998 0.102 0.042 0.203 0.409 0.133 0.048 0.306 0.591 0.087 0.026 0.198 0.409 0.152 0.069 0.314 0.591
		 1999 0.100 0.048 0.197 0.404 0.130 0.044 0.309 0.596 0.087 0.030 0.193 0.404 0.152 0.067 0.319 0.596
		 2000 0.104 0.044 0.207 0.452 0.127 0.047 0.288 0.548 0.091 0.026 0.202 0.452 0.146 0.066 0.297 0.548
		 2001 0.104 0.044 0.196 0.376 0.149 0.055 0.331 0.624 0.092 0.028 0.193 0.376 0.170 0.079 0.340 0.624
		 2002 0.106 0.044 0.218 0.471 0.107 0.039 0.263 0.529 0.091 0.027 0.213 0.471 0.123 0.055 0.274 0.529
		 2003 0.123 0.051 0.240 0.499 0.105 0.040 0.245 0.501 0.108 0.032 0.235 0.499 0.121 0.058 0.253 0.501
		 2004 0.102 0.041 0.200 0.419 0.112 0.035 0.277 0.581 0.089 0.026 0.195 0.419 0.132 0.055 0.285 0.581
		 2005 0.103 0.049 0.200 0.442 0.116 0.039 0.275 0.558 0.091 0.029 0.194 0.442 0.134 0.058 0.281 0.558
		 2006 0.111 0.050 0.209 0.398 0.125 0.047 0.290 0.602 0.097 0.030 0.203 0.398 0.143 0.067 0.300 0.602
		 2007 0.094 0.042 0.181 0.376 0.134 0.044 0.309 0.624 0.081 0.024 0.176 0.376 0.153 0.067 0.320 0.624
		 2008 0.107 0.050 0.212 0.465 0.120 0.044 0.267 0.535 0.095 0.031 0.208 0.465 0.134 0.061 0.277 0.535
		 2009 0.126 0.055 0.237 0.436 0.118 0.040 0.270 0.564 0.112 0.038 0.229 0.436 0.135 0.059 0.278 0.564
		 2010 0.112 0.045 0.223 0.471 0.107 0.039 0.245 0.529 0.095 0.028 0.216 0.471 0.123 0.058 0.254 0.529
	end
	gen ages="all"
	tempfile t1
	save `t1'
restore

* tabstat *ESI*_gYe [aw=weight], by(year) f(%4.3f)

clear
preserve
	input year ESI25_gYe ESI50_gYe ESI10_gYe ESI00_gYe upESI25_gYe upESI50_gYe upESI10_gYe upESI00_gYe gESI25_gYe gESI50_gYe gESI10_gYe gESI00_gYe upgESI25_gYe upgESI50_gYe upgESI10_gYe upgESI00_gYe
		1994 0.173 0.102 0.275 0.478 0.157 0.074 0.310 0.522 0.158 0.075 0.271 0.478 0.176 0.096 0.320 0.522
		1995 0.161 0.083 0.264 0.482 0.175 0.081 0.316 0.518 0.144 0.060 0.258 0.482 0.191 0.107 0.324 0.518
		1996 0.157 0.081 0.256 0.474 0.167 0.076 0.311 0.526 0.138 0.056 0.251 0.474 0.185 0.104 0.318 0.526
		1997 0.155 0.085 0.260 0.474 0.174 0.084 0.324 0.526 0.142 0.061 0.255 0.474 0.194 0.110 0.331 0.526
		1998 0.145 0.078 0.237 0.420 0.195 0.091 0.357 0.580 0.132 0.057 0.233 0.420 0.215 0.124 0.363 0.580
		1999 0.135 0.073 0.226 0.426 0.184 0.085 0.349 0.574 0.123 0.053 0.222 0.426 0.208 0.117 0.356 0.574
		2000 0.141 0.071 0.239 0.471 0.173 0.081 0.333 0.529 0.125 0.049 0.234 0.471 0.193 0.105 0.341 0.529
		2001 0.158 0.085 0.263 0.473 0.178 0.083 0.328 0.527 0.143 0.060 0.259 0.473 0.197 0.109 0.335 0.527
		2002 0.160 0.084 0.264 0.485 0.152 0.064 0.296 0.515 0.142 0.056 0.259 0.485 0.170 0.086 0.303 0.515
		2003 0.170 0.089 0.277 0.505 0.147 0.072 0.286 0.495 0.153 0.062 0.271 0.505 0.167 0.093 0.293 0.495
		2004 0.141 0.072 0.240 0.449 0.163 0.069 0.315 0.551 0.126 0.051 0.235 0.449 0.182 0.092 0.322 0.551
		2005 0.147 0.079 0.242 0.467 0.160 0.071 0.308 0.533 0.133 0.055 0.234 0.467 0.176 0.093 0.314 0.533
		2006 0.155 0.086 0.251 0.448 0.172 0.077 0.322 0.552 0.140 0.061 0.245 0.448 0.189 0.104 0.329 0.552
		2007 0.134 0.071 0.225 0.425 0.177 0.081 0.346 0.575 0.121 0.048 0.220 0.425 0.197 0.106 0.353 0.575
		2008 0.153 0.083 0.263 0.473 0.163 0.075 0.307 0.527 0.141 0.056 0.255 0.473 0.182 0.100 0.314 0.527
		2009 0.192 0.108 0.311 0.520 0.147 0.064 0.278 0.480 0.178 0.077 0.305 0.520 0.162 0.088 0.284 0.480
		2010 0.159 0.085 0.263 0.494 0.159 0.075 0.291 0.506 0.141 0.057 0.256 0.494 0.175 0.098 0.297 0.506
	end
	gen ages="all"
	tempfile t2
	save `t2'
restore

* tabstat RR??_?  [aw=weight], by(year) f(%4.3f)
clear
preserve
	input year RR00_1    RR00_2    RR00_3    RR00_4    RR10_1    RR10_2    RR10_3    RR10_4    RR25_1    RR25_2    RR25_3    RR25_4    RR50_1    RR50_2    RR50_3    RR50_4
		1994 0.447 0.075 0.092 0.386 0.685 0.040 0.073 0.201 0.808 0.019 0.061 0.112 0.891 0.008 0.044 0.057
		1995 0.434 0.084 0.075 0.407 0.700 0.037 0.068 0.196 0.824 0.016 0.063 0.098 0.912 0.006 0.040 0.042
		1996 0.463 0.063 0.084 0.390 0.709 0.035 0.075 0.181 0.829 0.014 0.068 0.089 0.915 0.004 0.044 0.037
		1997 0.469 0.058 0.085 0.388 0.709 0.031 0.079 0.181 0.833 0.013 0.060 0.095 0.910 0.005 0.047 0.038
		1998 0.517 0.062 0.078 0.342 0.734 0.030 0.065 0.171 0.844 0.012 0.059 0.086 0.920 0.002 0.043 0.035
		1999 0.512 0.062 0.093 0.332 0.741 0.034 0.065 0.160 0.851 0.015 0.054 0.080 0.920 0.007 0.037 0.036
		2000 0.474 0.055 0.079 0.392 0.727 0.034 0.065 0.174 0.845 0.014 0.054 0.087 0.926 0.004 0.034 0.036
		2001 0.496 0.031 0.127 0.346 0.719 0.018 0.083 0.179 0.835 0.007 0.061 0.096 0.913 0.002 0.044 0.040
		2002 0.465 0.051 0.075 0.410 0.711 0.025 0.073 0.191 0.831 0.008 0.065 0.095 0.913 0.003 0.045 0.038
		2003 0.435 0.060 0.067 0.438 0.693 0.031 0.067 0.209 0.818 0.012 0.061 0.108 0.908 0.004 0.043 0.045
		2004 0.499 0.053 0.084 0.364 0.735 0.026 0.065 0.175 0.846 0.014 0.054 0.087 0.924 0.004 0.037 0.035
		2005 0.484 0.048 0.083 0.385 0.735 0.024 0.064 0.177 0.844 0.009 0.054 0.093 0.918 0.003 0.035 0.044
		2006 0.509 0.043 0.092 0.356 0.727 0.022 0.064 0.187 0.836 0.009 0.055 0.101 0.911 0.003 0.042 0.045
		2007 0.528 0.048 0.098 0.326 0.754 0.022 0.066 0.158 0.857 0.009 0.052 0.082 0.925 0.004 0.036 0.035
		2008 0.476 0.051 0.065 0.408 0.711 0.026 0.075 0.187 0.838 0.009 0.056 0.097 0.912 0.005 0.040 0.043
		2009 0.446 0.035 0.113 0.406 0.669 0.020 0.094 0.217 0.799 0.008 0.077 0.115 0.889 0.002 0.058 0.050
		2010 0.451 0.055 0.082 0.413 0.710 0.028 0.068 0.195 0.829 0.012 0.062 0.096 0.912 0.004 0.046 0.038
	end
	gen ages="all"
	tempfile t3
	save `t3'
restore

* tabstat RR??_g? [aw=weight], by(year) f(%4.3f)
clear
preserve
	input year RR00_g1   RR00_g2   RR00_g3   RR00_g4   RR10_g1   RR10_g2   RR10_g3   RR10_g4   RR25_g1   RR25_g2   RR25_g3   RR25_g4   RR50_g1   RR50_g2   RR50_g3   RR50_g4
		1994 0.447 0.075 0.092 0.386 0.690 0.040 0.075 0.195 0.827 0.015 0.058 0.099 0.921 0.005 0.041 0.034
		1995 0.434 0.084 0.075 0.407 0.706 0.036 0.068 0.190 0.844 0.012 0.058 0.086 0.936 0.004 0.035 0.025
		1996 0.463 0.063 0.084 0.390 0.715 0.034 0.074 0.177 0.851 0.011 0.063 0.075 0.942 0.002 0.032 0.024
		1997 0.469 0.058 0.085 0.388 0.716 0.029 0.079 0.176 0.848 0.010 0.059 0.082 0.937 0.002 0.036 0.024
		1998 0.517 0.062 0.078 0.342 0.740 0.028 0.065 0.167 0.859 0.009 0.059 0.073 0.942 0.001 0.036 0.021
		1999 0.512 0.062 0.093 0.332 0.745 0.033 0.065 0.157 0.866 0.011 0.051 0.072 0.943 0.004 0.031 0.022
		2000 0.474 0.055 0.079 0.392 0.734 0.032 0.064 0.170 0.863 0.013 0.049 0.075 0.949 0.002 0.028 0.021
		2001 0.496 0.031 0.127 0.346 0.724 0.018 0.082 0.176 0.851 0.006 0.058 0.085 0.939 0.001 0.034 0.026
		2002 0.465 0.051 0.075 0.410 0.716 0.025 0.072 0.187 0.850 0.008 0.060 0.081 0.942 0.002 0.033 0.023
		2003 0.435 0.060 0.067 0.438 0.699 0.030 0.067 0.204 0.836 0.011 0.059 0.094 0.937 0.002 0.034 0.028
		2004 0.499 0.053 0.084 0.364 0.740 0.025 0.065 0.171 0.862 0.012 0.051 0.075 0.947 0.002 0.029 0.023
		2005 0.484 0.048 0.083 0.385 0.743 0.023 0.063 0.171 0.858 0.009 0.052 0.080 0.942 0.003 0.031 0.024
		2006 0.509 0.043 0.092 0.356 0.732 0.022 0.064 0.181 0.853 0.007 0.052 0.088 0.937 0.002 0.034 0.027
		2007 0.528 0.048 0.098 0.326 0.759 0.021 0.067 0.153 0.871 0.008 0.050 0.070 0.950 0.002 0.028 0.020
		2008 0.476 0.051 0.065 0.408 0.720 0.025 0.071 0.183 0.851 0.008 0.054 0.087 0.942 0.002 0.029 0.027
		2009 0.446 0.035 0.113 0.406 0.676 0.019 0.094 0.210 0.814 0.007 0.076 0.102 0.921 0.002 0.044 0.033
		2010 0.451 0.055 0.082 0.413 0.719 0.025 0.067 0.189 0.849 0.010 0.060 0.081 0.941 0.002 0.034 0.024
	end
	gen ages="all"
	tempfile t4
	save `t4'
restore
 
* tabstat d??_?Ye    [aw=weight], by(year) f(%4.3f) s(p50)
clear
preserve
	input year d00_nYe   d00_gYe   d10_nYe   d10_gYe   d25_nYe   d25_gYe   d50_nYe   d50_gYe
		1994 -0.111 -0.139 -0.279 -0.344 -0.519 -0.599 -0.712 -0.869
		1995 -0.093 -0.128 -0.250 -0.327 -0.448 -0.507 -0.767 -0.871
		1996 -0.093 -0.121 -0.246 -0.310 -0.450 -0.514 -0.787 -0.827
		1997 -0.089 -0.125 -0.261 -0.315 -0.452 -0.535 -0.763 -0.835
		1998 -0.099 -0.132 -0.250 -0.323 -0.435 -0.529 -0.773 -0.850
		1999 -0.094 -0.117 -0.252 -0.320 -0.483 -0.548 -0.760 -0.832
		2000 -0.084 -0.104 -0.252 -0.300 -0.444 -0.504 -0.737 -0.826
		2001 -0.110 -0.130 -0.268 -0.314 -0.454 -0.531 -0.731 -0.838
		2002 -0.086 -0.121 -0.241 -0.308 -0.441 -0.523 -0.752 -0.796
		2003 -0.095 -0.124 -0.256 -0.327 -0.434 -0.518 -0.761 -0.860
		2004 -0.092 -0.117 -0.253 -0.301 -0.432 -0.511 -0.755 -0.830
		2005 -0.083 -0.106 -0.259 -0.326 -0.476 -0.529 -0.726 -0.850
		2006 -0.109 -0.131 -0.266 -0.336 -0.468 -0.560 -0.734 -0.829
		2007 -0.093 -0.114 -0.259 -0.311 -0.457 -0.523 -0.713 -0.807
		2008 -0.084 -0.126 -0.250 -0.313 -0.475 -0.531 -0.732 -0.795
		2009 -0.117 -0.148 -0.269 -0.343 -0.454 -0.563 -0.765 -0.853
		2010 -0.091 -0.116 -0.251 -0.312 -0.434 -0.520 -0.780 -0.809
	end
	gen ages="all"
	tempfile t5
	save `t5'
restore

* tabstat d1_??_?Ye  [aw=weight], by(year) f(%4.3f) s(p50)
clear
preserve
	input year d1_00_nYe d1_00_gYe d1_10_nYe d1_10_gYe d1_25_nYe d1_25_gYe d1_50_nYe d1_50_gYe
		1994 -0.111 -0.139 -0.285 -0.356 -0.554 -0.639 -0.876 -1.039
		1995 -0.093 -0.128 -0.255 -0.332 -0.492 -0.557 -0.966 -1.045
		1996 -0.093 -0.121 -0.252 -0.318 -0.492 -0.569 -0.931 -1.005
		1997 -0.089 -0.125 -0.267 -0.321 -0.491 -0.581 -0.941 -1.019
		1998 -0.099 -0.132 -0.257 -0.328 -0.487 -0.578 -0.955 -1.029
		1999 -0.094 -0.117 -0.259 -0.326 -0.529 -0.595 -0.959 -1.027
		2000 -0.084 -0.104 -0.257 -0.307 -0.485 -0.559 -0.913 -1.018
		2001 -0.110 -0.130 -0.274 -0.321 -0.486 -0.568 -0.869 -1.013
		2002 -0.086 -0.121 -0.247 -0.316 -0.483 -0.574 -0.963 -1.007
		2003 -0.095 -0.124 -0.262 -0.333 -0.482 -0.553 -0.917 -1.050
		2004 -0.092 -0.117 -0.262 -0.307 -0.477 -0.571 -0.887 -0.991
		2005 -0.083 -0.106 -0.270 -0.337 -0.519 -0.579 -0.929 -1.016
		2006 -0.109 -0.131 -0.275 -0.345 -0.506 -0.598 -0.930 -1.014
		2007 -0.093 -0.114 -0.264 -0.318 -0.514 -0.578 -0.915 -0.979
		2008 -0.084 -0.126 -0.257 -0.326 -0.508 -0.562 -0.907 -1.055
		2009 -0.117 -0.148 -0.279 -0.352 -0.492 -0.601 -0.892 -1.021
		2010 -0.091 -0.116 -0.259 -0.320 -0.488 -0.572 -0.975 -1.033
	end
	gen ages="all"
	tempfile t6
	save `t6'
restore


**************
* ESIs, ages 25-60
**************

* tabstat *ESI*_nYe [aw=weight] if age>=25 & age<=60, by(year) f(%4.3f)

clear
preserve
	input year ESI25_nYe ESI50_nYe ESI10_nYe ESI00_nYe upESI25_nYe upESI50_nYe upESI10_nYe upESI00_nYe gESI25_nYe gESI50_nYe gESI10_nYe gESI00_nYe upgESI25_nYe upgESI50_nYe upgESI10_nYe upgESI00_nYe 
		1994 0.146 0.072 0.257 0.471 0.109 0.038 0.267 0.529 0.128 0.044 0.250 0.471 0.126 0.055 0.276 0.529
		1995 0.128 0.053 0.254 0.507 0.123 0.047 0.264 0.493 0.111 0.032 0.248 0.507 0.139 0.068 0.272 0.493
		1996 0.116 0.047 0.233 0.457 0.118 0.043 0.270 0.543 0.097 0.030 0.228 0.457 0.134 0.061 0.277 0.543
		1997 0.117 0.049 0.225 0.451 0.121 0.044 0.281 0.549 0.101 0.030 0.219 0.451 0.140 0.065 0.289 0.549
		1998 0.108 0.045 0.211 0.408 0.135 0.049 0.310 0.592 0.092 0.028 0.206 0.408 0.154 0.070 0.319 0.592
		1999 0.107 0.050 0.208 0.407 0.128 0.043 0.314 0.593 0.093 0.032 0.203 0.407 0.153 0.065 0.325 0.593
		2000 0.109 0.044 0.217 0.448 0.131 0.046 0.296 0.552 0.095 0.026 0.212 0.448 0.150 0.066 0.306 0.552
		2001 0.111 0.044 0.208 0.387 0.152 0.056 0.338 0.613 0.098 0.027 0.204 0.387 0.174 0.078 0.347 0.613
		2002 0.109 0.044 0.227 0.460 0.110 0.038 0.271 0.540 0.093 0.027 0.222 0.460 0.127 0.055 0.281 0.540
		2003 0.130 0.055 0.250 0.498 0.104 0.039 0.248 0.502 0.115 0.033 0.245 0.498 0.120 0.056 0.255 0.502
		2004 0.109 0.043 0.213 0.428 0.115 0.037 0.283 0.572 0.096 0.027 0.208 0.428 0.136 0.056 0.291 0.572
		2005 0.110 0.052 0.211 0.434 0.123 0.041 0.286 0.566 0.097 0.031 0.204 0.434 0.142 0.062 0.291 0.566
		2006 0.120 0.052 0.223 0.420 0.125 0.047 0.294 0.580 0.105 0.031 0.217 0.420 0.145 0.066 0.304 0.580
		2007 0.103 0.045 0.194 0.389 0.138 0.045 0.316 0.611 0.089 0.026 0.188 0.389 0.158 0.070 0.328 0.611
		2008 0.114 0.053 0.227 0.457 0.120 0.043 0.274 0.543 0.102 0.031 0.223 0.457 0.136 0.060 0.284 0.543
		2009 0.136 0.058 0.254 0.450 0.121 0.042 0.277 0.550 0.122 0.040 0.245 0.450 0.139 0.060 0.285 0.550
		2010 0.118 0.048 0.234 0.461 0.113 0.041 0.256 0.539 0.100 0.028 0.227 0.461 0.129 0.060 0.265 0.539
	end
	gen ages="25-60"
	tempfile T1
	save `T1'
restore

* tabstat *ESI*_gYe [aw=weight] if age>=25 & age<=60, by(year) f(%4.3f)

clear
preserve
	input year ESI25_gYe ESI50_gYe ESI10_gYe ESI00_gYe upESI25_gYe upESI50_gYe upESI10_gYe upESI00_gYe gESI25_gYe gESI50_gYe gESI10_gYe gESI00_gYe upgESI25_gYe upgESI50_gYe upgESI10_gYe upgESI00_gYe
		1994 0.176 0.102 0.271 0.468 0.159 0.071 0.318 0.532 0.159 0.074 0.267 0.468 0.178 0.095 0.327 0.532
		1995 0.162 0.080 0.270 0.486 0.171 0.078 0.314 0.514 0.146 0.056 0.265 0.486 0.187 0.103 0.321 0.514
		1996 0.159 0.082 0.260 0.468 0.168 0.074 0.315 0.532 0.141 0.055 0.254 0.468 0.185 0.104 0.324 0.532
		1997 0.147 0.077 0.251 0.457 0.174 0.080 0.328 0.543 0.134 0.053 0.247 0.457 0.193 0.107 0.335 0.543
		1998 0.141 0.073 0.233 0.410 0.188 0.083 0.358 0.590 0.127 0.051 0.228 0.410 0.208 0.115 0.365 0.590
		1999 0.134 0.071 0.227 0.427 0.175 0.078 0.342 0.573 0.121 0.048 0.224 0.427 0.199 0.108 0.349 0.573
		2000 0.139 0.065 0.236 0.462 0.170 0.076 0.338 0.538 0.123 0.043 0.230 0.462 0.190 0.099 0.346 0.538
		2001 0.157 0.080 0.265 0.470 0.176 0.077 0.327 0.530 0.142 0.054 0.260 0.470 0.195 0.104 0.334 0.530
		2002 0.152 0.075 0.253 0.467 0.151 0.059 0.302 0.533 0.131 0.048 0.249 0.467 0.169 0.081 0.310 0.533
		2003 0.165 0.084 0.274 0.497 0.142 0.066 0.287 0.503 0.149 0.057 0.268 0.497 0.162 0.087 0.294 0.503
		2004 0.141 0.068 0.243 0.441 0.162 0.065 0.321 0.559 0.125 0.047 0.239 0.441 0.182 0.088 0.328 0.559
		2005 0.146 0.077 0.239 0.453 0.162 0.071 0.316 0.547 0.131 0.052 0.232 0.453 0.179 0.094 0.322 0.547
		2006 0.156 0.083 0.255 0.445 0.168 0.073 0.321 0.555 0.141 0.056 0.248 0.445 0.186 0.098 0.328 0.555
		2007 0.138 0.072 0.229 0.423 0.176 0.078 0.345 0.576 0.124 0.046 0.224 0.423 0.196 0.103 0.352 0.576
		2008 0.152 0.078 0.256 0.454 0.161 0.071 0.311 0.546 0.138 0.052 0.249 0.454 0.180 0.096 0.318 0.546
		2009 0.193 0.106 0.313 0.511 0.146 0.060 0.285 0.489 0.178 0.074 0.306 0.511 0.162 0.085 0.291 0.489
		2010 0.156 0.080 0.259 0.475 0.160 0.072 0.297 0.525 0.137 0.053 0.252 0.475 0.176 0.096 0.303 0.525
	end
	gen ages="25-60"
	tempfile T2
	save `T2'
restore

* tabstat RR??_? [aw=weight] if age>=25 & age<=60, by(year) f(%4.3f)

clear
preserve
	input year RR00_1 RR00_2 RR00_3 RR00_4 RR10_1 RR10_2 RR10_3 RR10_4 RR25_1 RR25_2 RR25_3 RR25_4 RR50_1 RR50_2 RR50_3 RR50_4
		1994 0.456 0.077 0.075 0.393 0.687 0.042 0.059 0.211 0.804 0.021 0.055 0.120 0.889 0.009 0.042 0.059
		1995 0.435 0.079 0.066 0.420 0.692 0.038 0.061 0.209 0.820 0.019 0.059 0.103 0.913 0.007 0.039 0.041
		1996 0.472 0.059 0.075 0.394 0.703 0.037 0.068 0.192 0.827 0.014 0.062 0.097 0.914 0.004 0.042 0.040
		1997 0.484 0.059 0.073 0.384 0.717 0.032 0.061 0.190 0.840 0.014 0.048 0.098 0.918 0.006 0.038 0.038
		1998 0.530 0.060 0.069 0.341 0.738 0.030 0.054 0.178 0.847 0.012 0.050 0.091 0.925 0.003 0.036 0.036
		1999 0.514 0.059 0.088 0.338 0.736 0.037 0.059 0.167 0.850 0.017 0.047 0.085 0.923 0.007 0.032 0.038
		2000 0.487 0.052 0.070 0.391 0.730 0.034 0.055 0.181 0.847 0.015 0.049 0.089 0.931 0.004 0.030 0.036
		2001 0.500 0.030 0.113 0.356 0.717 0.018 0.074 0.190 0.836 0.008 0.055 0.102 0.918 0.002 0.039 0.040
		2002 0.483 0.050 0.066 0.400 0.720 0.026 0.056 0.198 0.839 0.009 0.055 0.097 0.922 0.003 0.037 0.038
		2003 0.444 0.059 0.062 0.435 0.695 0.032 0.055 0.219 0.821 0.014 0.051 0.114 0.911 0.005 0.036 0.048
		2004 0.506 0.053 0.070 0.372 0.731 0.026 0.056 0.187 0.844 0.015 0.049 0.092 0.928 0.005 0.032 0.036
		2005 0.500 0.047 0.073 0.380 0.737 0.024 0.053 0.186 0.845 0.009 0.047 0.099 0.920 0.003 0.032 0.046
		2006 0.511 0.043 0.071 0.374 0.722 0.022 0.056 0.199 0.834 0.010 0.049 0.107 0.913 0.003 0.037 0.046
		2007 0.532 0.045 0.083 0.340 0.750 0.021 0.059 0.170 0.852 0.010 0.049 0.089 0.924 0.004 0.035 0.036
		2008 0.496 0.051 0.053 0.401 0.715 0.030 0.060 0.196 0.839 0.010 0.049 0.102 0.917 0.005 0.032 0.045
		2009 0.455 0.035 0.093 0.417 0.667 0.020 0.080 0.232 0.798 0.009 0.069 0.124 0.892 0.003 0.053 0.052
		2010 0.468 0.057 0.071 0.404 0.712 0.030 0.055 0.203 0.831 0.013 0.054 0.102 0.916 0.004 0.039 0.040
	end
	gen ages="25-60"
	tempfile T3
	save `T3'
restore

* tabstat RR??_g? [aw=weight] if age>=25 & age<=60, by(year) f(%4.3f)

clear
preserve
	input year RR00_g1 RR00_g2 RR00_g3 RR00_g4 RR10_g1 RR10_g2 RR10_g3 RR10_g4 RR25_g1 RR25_g2 RR25_g3 RR25_g4 RR50_g1 RR50_g2 RR50_g3 RR50_g4
		1994 0.456 0.077 0.075 0.393 0.692 0.042 0.063 0.204 0.825 0.017 0.052 0.106 0.921 0.005 0.038 0.036
		1995 0.435 0.079 0.066 0.420 0.698 0.037 0.061 0.204 0.840 0.015 0.056 0.089 0.939 0.004 0.033 0.024
		1996 0.472 0.059 0.075 0.394 0.709 0.037 0.067 0.188 0.848 0.011 0.060 0.082 0.943 0.002 0.031 0.024
		1997 0.484 0.059 0.073 0.384 0.723 0.030 0.061 0.185 0.855 0.011 0.049 0.085 0.945 0.002 0.029 0.024
		1998 0.530 0.060 0.069 0.341 0.743 0.029 0.054 0.174 0.863 0.010 0.050 0.076 0.948 0.002 0.029 0.021
		1999 0.514 0.059 0.088 0.338 0.741 0.036 0.059 0.164 0.867 0.012 0.044 0.076 0.948 0.004 0.025 0.023
		2000 0.487 0.052 0.070 0.391 0.737 0.033 0.053 0.177 0.865 0.013 0.045 0.077 0.955 0.002 0.023 0.020
		2001 0.500 0.030 0.113 0.356 0.722 0.018 0.073 0.186 0.852 0.006 0.052 0.089 0.945 0.001 0.029 0.025
		2002 0.483 0.050 0.066 0.400 0.725 0.026 0.056 0.193 0.861 0.009 0.049 0.082 0.950 0.002 0.026 0.023
		2003 0.444 0.059 0.062 0.435 0.700 0.032 0.055 0.213 0.839 0.012 0.049 0.100 0.941 0.002 0.028 0.029
		2004 0.506 0.053 0.070 0.372 0.736 0.026 0.057 0.182 0.862 0.014 0.045 0.080 0.950 0.002 0.025 0.023
		2005 0.500 0.047 0.073 0.380 0.746 0.023 0.051 0.181 0.859 0.010 0.047 0.084 0.945 0.003 0.027 0.025
		2006 0.511 0.043 0.071 0.374 0.728 0.023 0.056 0.192 0.851 0.008 0.047 0.094 0.942 0.002 0.030 0.026
		2007 0.532 0.045 0.083 0.340 0.755 0.021 0.059 0.165 0.867 0.009 0.048 0.076 0.951 0.003 0.025 0.021
		2008 0.496 0.051 0.053 0.401 0.723 0.028 0.056 0.193 0.853 0.009 0.046 0.092 0.945 0.002 0.026 0.026
		2009 0.455 0.035 0.093 0.417 0.675 0.019 0.081 0.225 0.814 0.008 0.067 0.111 0.924 0.002 0.039 0.034
		2010 0.468 0.057 0.071 0.404 0.720 0.027 0.054 0.198 0.851 0.012 0.053 0.084 0.945 0.002 0.030 0.023
	end
	gen ages="25-60"
	tempfile T4
	save `T4'
restore
 
* tabstat d??_?Ye [aw=weight] if age>=25 & age<=60, by(year) f(%4.3f) s(p50)
clear
preserve
	input year d00_nYe d00_gYe d10_nYe d10_gYe d25_nYe d25_gYe d50_nYe d50_gYe
		1994 -0.118 -0.145 -0.285 -0.357 -0.487 -0.583 -0.713 -0.869
		1995 -0.100 -0.133 -0.252 -0.320 -0.445 -0.491 -0.759 -0.847
		1996 -0.105 -0.128 -0.249 -0.316 -0.444 -0.514 -0.764 -0.799
		1997 -0.100 -0.125 -0.259 -0.307 -0.446 -0.509 -0.769 -0.817
		1998 -0.108 -0.136 -0.255 -0.312 -0.440 -0.513 -0.765 -0.836
		1999 -0.105 -0.119 -0.255 -0.309 -0.478 -0.529 -0.763 -0.783
		2000 -0.093 -0.105 -0.251 -0.299 -0.439 -0.473 -0.733 -0.784
		2001 -0.114 -0.133 -0.269 -0.309 -0.446 -0.510 -0.722 -0.801
		2002 -0.097 -0.121 -0.241 -0.297 -0.439 -0.491 -0.741 -0.771
		2003 -0.101 -0.125 -0.261 -0.322 -0.433 -0.506 -0.748 -0.837
		2004 -0.100 -0.125 -0.257 -0.293 -0.428 -0.485 -0.743 -0.803
		2005 -0.098 -0.111 -0.262 -0.324 -0.474 -0.528 -0.726 -0.811
		2006 -0.113 -0.136 -0.272 -0.331 -0.458 -0.532 -0.734 -0.782
		2007 -0.100 -0.120 -0.261 -0.313 -0.452 -0.517 -0.716 -0.795
		2008 -0.099 -0.131 -0.250 -0.312 -0.478 -0.507 -0.725 -0.765
		2009 -0.126 -0.156 -0.272 -0.336 -0.440 -0.545 -0.767 -0.824
		2010 -0.102 -0.121 -0.252 -0.304 -0.429 -0.510 -0.745 -0.789
	end
	gen ages="25-60"
	tempfile T5
	save `T5'
restore

* tabstat d1_??_?Ye [aw=weight] if age>=25 & age<=60, by(year) f(%4.3f) s(p50)

clear
preserve
	input year d1_00_nYe d1_00_gYe d1_10_nYe d1_10_gYe d1_25_nYe d1_25_gYe d1_50_nYe d1_50_gYe
		1994 -0.118 -0.145 -0.294 -0.364 -0.549 -0.633 -0.872 -1.051
		1995 -0.100 -0.133 -0.259 -0.327 -0.483 -0.542 -0.989 -1.035
		1996 -0.105 -0.128 -0.255 -0.323 -0.492 -0.563 -0.931 -0.981
		1997 -0.100 -0.125 -0.264 -0.312 -0.485 -0.548 -0.941 -1.007
		1998 -0.108 -0.136 -0.260 -0.317 -0.487 -0.551 -0.963 -1.034
		1999 -0.105 -0.119 -0.263 -0.314 -0.528 -0.572 -0.928 -1.004
		2000 -0.093 -0.105 -0.258 -0.306 -0.473 -0.526 -0.901 -0.998
		2001 -0.114 -0.133 -0.274 -0.314 -0.477 -0.543 -0.856 -0.996
		2002 -0.097 -0.121 -0.246 -0.301 -0.481 -0.561 -0.963 -0.982
		2003 -0.101 -0.125 -0.268 -0.327 -0.483 -0.538 -0.910 -1.039
		2004 -0.100 -0.125 -0.264 -0.299 -0.466 -0.549 -0.865 -0.976
		2005 -0.098 -0.111 -0.274 -0.338 -0.519 -0.573 -0.904 -1.001
		2006 -0.113 -0.136 -0.278 -0.338 -0.499 -0.576 -0.957 -0.996
		2007 -0.100 -0.120 -0.270 -0.319 -0.506 -0.557 -0.908 -0.963
		2008 -0.099 -0.131 -0.256 -0.326 -0.508 -0.542 -0.907 -0.998
		2009 -0.126 -0.156 -0.282 -0.345 -0.481 -0.576 -0.895 -1.001
		2010 -0.102 -0.121 -0.259 -0.311 -0.488 -0.561 -0.976 -1.003
	end
	gen ages="25-60"
	tempfile T6
	save `T6'
restore

clear
set obs 1
gen year=-99
gen ages="junk"
forvalues j=1/6 {
	cap drop _merge
	merge 1:1 year ages using `t`j''
}
tempfile t
save `t'

clear
set obs 1
gen year=-99
gen ages="junk"
forvalues j=1/6 {
	cap drop _merge
	merge 1:1 year ages using `T`j''
}
append using `t'
drop if year==-99
isid year ages

label var ESI25_nYe "25+% arc drop, HH disposable income"
label var ESI50_nYe "50+% arc drop, HH disposable income"
label var ESI10_nYe "10+% arc drop, HH disposable income"
label var ESI00_nYe "00+% arc drop, HH disposable income"
label var upESI25_nYe "25+% arc gain, HH disposable income"
label var upESI50_nYe "50+% arc gain, HH disposable income"
label var upESI10_nYe "10+% arc gain, HH disposable income"
label var upESI00_nYe "00+% arc gain, HH disposable income"
label var gESI25_nYe "25+% drop, HH disposable income"
label var gESI50_nYe "50+% drop, HH disposable income"
label var gESI10_nYe "10+% drop, HH disposable income"
label var gESI00_nYe "00+% drop, HH disposable income"
label var upgESI25_nYe "25+% gain, HH disposable income"
label var upgESI50_nYe "50+% gain, HH disposable income"
label var upgESI10_nYe "10+% gain, HH disposable income"
label var upgESI00_nYe "00+% gain, HH disposable income"
label var ESI25_gYe "25+% arc drop, HH market income"
label var ESI50_gYe "50+% arc drop, HH market income"
label var ESI10_gYe "10+% arc drop, HH market income"
label var ESI00_gYe "00+% arc drop, HH market income"
label var upESI25_gYe "25+% arc gain, HH market income"
label var upESI50_gYe "50+% arc gain, HH market income"
label var upESI10_gYe "10+% arc gain, HH market income"
label var upESI00_gYe "00+% arc gain, HH market income"
label var gESI25_gYe "25+% drop, HH market income"
label var gESI50_gYe "50+% drop, HH market income"
label var gESI10_gYe "10+% drop, HH market income"
label var gESI00_gYe "00+% drop, HH market income"
label var upgESI25_gYe "25+% gain, HH market income"
label var upgESI50_gYe "50+% gain, HH market income"
label var upgESI10_gYe "10+% gain, HH market income"
label var upgESI00_gYe "00+% gain, HH market income"
label var RR00_1 "MI+DI drops (ESI00_gYe=0 + ESI00_nYe==0)"
label var RR00_2 "MI+DI drops (ESI00_gYe=0 + ESI00_nYe==1)"
label var RR00_3 "MI+DI drops (ESI00_gYe=1 + ESI00_nYe==0)"
label var RR00_4 "MI+DI drops (ESI00_gYe=1 + ESI00_nYe==1)"
label var RR10_1 "MI+DI drops (ESI10_gYe=0 + ESI10_nYe==0)"
label var RR10_2 "MI+DI drops (ESI10_gYe=0 + ESI10_nYe==1)"
label var RR10_3 "MI+DI drops (ESI10_gYe=1 + ESI10_nYe==0)"
label var RR10_4 "MI+DI drops (ESI10_gYe=1 + ESI10_nYe==1)"
label var RR25_1 "MI+DI drops (ESI25_gYe=0 + ESI25_nYe==0)"
label var RR25_2 "MI+DI drops (ESI25_gYe=0 + ESI25_nYe==1)"
label var RR25_3 "MI+DI drops (ESI25_gYe=1 + ESI25_nYe==0)"
label var RR25_4 "MI+DI drops (ESI25_gYe=1 + ESI25_nYe==1)"
label var RR50_1 "MI+DI drops (ESI50_gYe=0 + ESI50_nYe==0)"
label var RR50_2 "MI+DI drops (ESI50_gYe=0 + ESI50_nYe==1)"
label var RR50_3 "MI+DI drops (ESI50_gYe=1 + ESI50_nYe==0)"
label var RR50_4 "MI+DI drops (ESI50_gYe=1 + ESI50_nYe==1)"
label var RR00_g1 "MI+DI drops (gESI00_gYe=0 + gESI00_nYe==0)"
label var RR00_g2 "MI+DI drops (gESI00_gYe=0 + gESI00_nYe==1)"
label var RR00_g3 "MI+DI drops (gESI00_gYe=1 + gESI00_nYe==0)"
label var RR00_g4 "MI+DI drops (gESI00_gYe=1 + gESI00_nYe==1)"
label var RR10_g1 "MI+DI drops (gESI10_gYe=0 + gESI10_nYe==0)"
label var RR10_g2 "MI+DI drops (gESI10_gYe=0 + gESI10_nYe==1)"
label var RR10_g3 "MI+DI drops (gESI10_gYe=1 + gESI10_nYe==0)"
label var RR10_g4 "MI+DI drops (gESI10_gYe=1 + gESI10_nYe==1)"
label var RR25_g1 "MI+DI drops (gESI25_gYe=0 + gESI25_nYe==0)"
label var RR25_g2 "MI+DI drops (gESI25_gYe=0 + gESI25_nYe==1)"
label var RR25_g3 "MI+DI drops (gESI25_gYe=1 + gESI25_nYe==0)"
label var RR25_g4 "MI+DI drops (gESI25_gYe=1 + gESI25_nYe==1)"
label var RR50_g1 "MI+DI drops (gESI50_gYe=0 + gESI50_nYe==0)"
label var RR50_g2 "MI+DI drops (gESI50_gYe=0 + gESI50_nYe==1)"
label var RR50_g3 "MI+DI drops (gESI50_gYe=1 + gESI50_nYe==0)"
label var RR50_g4 "MI+DI drops (gESI50_gYe=1 + gESI50_nYe==1)"
label var d00_nYe "Median drop size if ESI00_nYe==1"
label var d00_gYe "Median drop size if ESI00_gYe==1"
label var d10_nYe "Median drop size if ESI10_nYe==1"
label var d10_gYe "Median drop size if ESI10_gYe==1"
label var d25_nYe "Median drop size if ESI25_nYe==1"
label var d25_gYe "Median drop size if ESI25_gYe==1"
label var d50_nYe "Median drop size if ESI50_nYe==1"
label var d50_gYe "Median drop size if ESI50_gYe==1"
label var d1_00_nYe "Median drop size if gESI00_nYe==1"
label var d1_00_gYe "Median drop size if gESI00_gYe==1"
label var d1_10_nYe "Median drop size if gESI10_nYe==1"
label var d1_10_gYe "Median drop size if gESI10_gYe==1"
label var d1_25_nYe "Median drop size if gESI25_nYe==1"
label var d1_25_gYe "Median drop size if gESI25_gYe==1"
label var d1_50_nYe "Median drop size if gESI50_nYe==1"
label var d1_50_gYe "Median drop size if gESI50_gYe==1"

* gini coefficients
* foreach v in nY gY nYe gYe {
*	di "** zeros excluded, ages 25-60 **"
*	ineqdec0 `v' [aw=weight] if age>=25 & age<=60, bygroup(year)
*	di "** zeros excluded, all ages **"
*	ineqdec0 `v' [aw=weight], bygroup(year)
*	replace `v'=1 if `v'>=0 & `v'<1 // for well-defined ginis. Do we want to do this??
*	di "** income bottom coded at 1, ages 25-60 **"
*	ineqdec0 `v' [aw=weight] if age>=25 & age<=60, bygroup(year)
*	di "** income bottom coded at 1, all ages **"
*	ineqdec0 `v' [aw=weight], bygroup(year)
* }

preserve
	clear
	input year	str7 ages gini_nY gini_gY gini_nYe gini_gYe
		1993 25-60 0.30632 0.37436 0.28295 0.35767
		1994 25-60 0.31073 0.37606 0.28717 0.36083
		1995 25-60 0.30976 0.37633 0.28338 0.35454
		1996 25-60 0.31875 0.38534 0.29144 0.36322
		1997 25-60 0.32731 0.39242 0.29591 0.36588
		1998 25-60 0.32937 0.39154 0.29518 0.36185
		1999 25-60 0.33077 0.39144 0.29273 0.35559
		2000 25-60 0.33629 0.39547 0.297 0.35945
		2001 25-60 0.33951 0.39991 0.29912 0.36386
		2002 25-60 0.33794 0.39908 0.29862 0.36278
		2003 25-60 0.33396 0.39429 0.29886 0.36334
		2004 25-60 0.34265 0.40347 0.30265 0.36688
		2005 25-60 0.33779 0.39528 0.30021 0.36365
		2006 25-60 0.33687 0.39635 0.30147 0.36713
		2007 25-60 0.3353 0.39753 0.29853 0.36723
		2008 25-60 0.33747 0.39704 0.29755 0.36213
		2009 25-60 0.33738 0.40021 0.30124 0.3704
		2010 25-60 0.33587 0.39839 0.29884 0.36709
		1993 all 0.3152 0.40868 0.2848 0.38855
		1994 all 0.31998 0.41041 0.2887 0.39047
		1995 all 0.31845 0.41042 0.28567 0.38599
		1996 all 0.32672 0.41861 0.29344 0.3942
		1997 all 0.33326 0.42311 0.29651 0.39573
		1998 all 0.33503 0.42403 0.29603 0.39389
		1999 all 0.33604 0.42217 0.29444 0.38855
		2000 all 0.34488 0.42944 0.29925 0.39203
		2001 all 0.34511 0.43087 0.30143 0.39623
		2002 all 0.34411 0.42923 0.30192 0.39513
		2003 all 0.34197 0.42764 0.30159 0.39527
		2004 all 0.34962 0.43618 0.30514 0.39914
		2005 all 0.34552 0.42862 0.30257 0.3951
		2006 all 0.34185 0.42783 0.30115 0.3964
		2007 all 0.34057 0.42771 0.30039 0.39768
		2008 all 0.3474 0.43241 0.30191 0.3957
		2009 all 0.34542 0.43467 0.30299 0.40213
		2010 all 0.34592 0.43634 0.30146 0.40128
	end
	isid year ages
	foreach v in gini_nY gini_gY gini_nYe gini_gYe {
		gen `v'_all=`v' if ages=="all"
		gen `v'_2560=`v' if ages=="25-60"
		drop `v'
	}
	label var gini_nYe_2560	"Gini if HH disposable income, ages 25-60"
	label var gini_nYe_all 	"Gini if HH disposable income, ages all"
	label var gini_gYe_2560	"Gini if HH market income, ages 25-60"
	label var gini_gYe_all 	"Gini if HH market income, ages all"
	drop gini_nY_all gini_nY_2560 gini_gY_all gini_gY_2560
	tempfile t
	save `t'
restore

cap drop _merge
merge 1:1 year ages using `t'
drop _merge

compress
gen dataset="CNEF"
gen source="SLID"
gen ccode=20
note: C:\Dropbox\Rehm\Data\CNEF\cnef_1970-2014\slid\slid_rehm_v9_log_for-stata.do
note: Created on `= c(current_date)'

saveold slid_1993-2010.dta, replace

